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Magnetic Fields: Reasons for Simulation and Methods Available 


Norman F. Ness 

, . . . . NASA x Goddard Space Flight Center 

Introduction 

Magnetic fields are among the most ubiquitous of physical 
phenomenon studied by man. They are important on an atomic 
level in which magnetic field strengths of thousands of gauss are 
employed in the study of atomic structure. In cosmic physics 
stellar magnetic fields of tens and hundreds of gauss are 
important factors controlling the physical phenomenon observed. 

In interstellar and in interplanetary space magnetic fields of 
less than 10 - ^ gauss are important in determining the motion of 
charged particles and in the dynamics of the interstellar and 
interplanetary medium. The major portion of this talk is to 
be directed towards a description of the present spatial 
environment of the Earth with respect to magnetic fields. 

Following the descriptive aspects of the magnetic field as observed 
by satellite and space probes, the quantitative representation 
of these fields will be given. Finally a description of the 
presently available systems to simulate the magnetic fields in 
space will be reviewed and the most recent large scale facility 
for such work discussed. 

The last section of the paper deals with engineering 

applications of our knowledge of the earth’s magnetic field in 

various satellite programs. Certain of these applications are 

quite unique and represent some of the reasons for necessitating 

simulation of magnetic fields as observed in space. The unit of 

magnetic field force as commonly used in geomagnetism is the 
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gamma (equal to 10 Gauss) . On this scale a current of one 
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milliampere flowing in an infinitely long wire creates a 
magnetic field whose strength is 20 gammas at a distance 
of one centimeter. 

There are a number of scientific reasons for mapping the 
magnetic field environment of the earth. These include: 

1. Investigation of secular changes and the present 

state so as to determine the sources of the geomagnetic 
field and the interplanetary magnetic field structure; 

2. Tb study and analyse charged particle motion which 
requires a knowledge of the geomagnetic and inter- 
planetary magnetic fields, and 

3. To determine the time variations of these fields 
and their correlation with auroral phenomenon, and 
solar activity. 

In order to successfully map the magnetic fields in space 
it is necessary to perform direct measurements of the magnetic 
field from satellites and space probes. Thus they represent 
portable laboratory platforms upon which instruments sensitive 
to magnetic fields are placed. One of the primary engineering 
reasons for simulating magnetic fields in space is to test and 
calibrate the magnetometer sensors which will be used to directly 
measure magnetic fields in space. A second objective is to 
allow a measurement of the magnetic properties of satellite and 
space probes to determine : 
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1. Their contaminating or noise effects on magnetic 
field measurements made by magnetometers on board 
the same spacecraft, 

2. To determine the dynamic interaction of the moving 
and possibly spinning spacecraft with magnetic fields, 
principally the earth’s field and 

3. To allow the testing of specific orientation and 
attitude control and spin producing subsystems which 
employ the geomagnetic field for magnetic torque 
interactions and/or directional information. 

As more is learned about the magnetic fields in space more 
sophisticated use in an engineering sense can be made of this 
valuable information. The principal techniques of simulating 
terrestrially the weak magnetic fields which exist in space employ 
large current carrying coil systems which generate steady 
magnetic fields canceling the earth’s field over a small volume. 
Results of Magnetic Field Measurements in Space 

A number of satellites carrying magnetometers have been 
launched by both this country and the USSR. A summary of those 
satellites measuring principally the earth’s magnetic field is 
shown in Figure 1. In general the life time of these satellites 
has been individually limited as shown in Figure 2, summarizing 
trajectory characteristics pertinent to mapping of the earth’s 
magnetic field. A very limited number of space probes have 
investigated the interplanetary magnetic field. In order to 
measure the properties of the interplanetary medium, undisturbed 
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by the presence of the Earth, it is necessary to perform 
measurements at a considerable distance from the surface of 
the earth. The interaction of the solar plasma resulting from 
the expanding solar corona requires that these measurements be 
performed at least 80,000 Km from the earth near the subsolar 
point and at a considerably larger distance away from this 
region. As shown in Figure 3 those satellites which have 
performed measurements of the interplanetary field were quite 
limited either in accuracy or in measuring only a limited 
characteristic of the field and not the complete vector magnetic 
f ield. 

A representative sample of the results of measurements in 
^>ace confirming the general characteristics of the dipolar 

< 

geomagnetic field is shown in Figure 4. The Explorer X satellite 
carried a rubidium vapor magnetometer and performed successful 
measurements of the earth’s magnetic field from 2 to 7 Rg (Earth 
radii) in March 1961. In Figure 4 are summarized both the direct 
measurements and the differences between the observations and 
the theoretical field values predicted by extrapolation using 
spherical harmonic analyses of the surface terrestrial field. 

It is seen that the deviations are quite small and represent less 
than 5 percent of the ambient field. However, these small 
deviations of the earth's magnetic field are extremely important 
in this region of space. It is precisely within this radial 
distance range that the trapped charged particle fluxes are 
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observed to reach their peak intensities. The deviations 

observed, as illustrated in Figure 4, are consistant with 

the magnetic effects of charged particle spiral and drift 

motion in ±he earth's magnetic field. At the present time 

detailed models of the particle belts are being developed as 

more refined instrumentation permits investigation of the 

spectral and pitch angle distributions of energetic particles 

4 

within the belts. 

A substantially different result in measuring the 
earth's magnetic field results when one considers Figure 5. 

5 

The magnetic field results as obtained on the IMP-1 satellite 

are shown near the subsolar point. It is seen that the earth’s 

magnetic field agrees reasonably well with that predicted by 

spherical harmonic analysis until a distance of 7 or 8 R e is 

reached. Beyond this point the observed field becomes increasingly 

stronger until a distance of 10.7 R 0 is reached. At this point 

the measured magnitude of the field is approximately twice that 

predicted by theoretical analysis. Subsequent to this point on 

the trajectory the field abruptly drops to low and fluctuating 

values. This characteristic abrupt decrease of the earth’s 

magnetic field following a gradual but significant increase in 

magnetic is identified as a boundary of the earth's magnetic 

field. It is due to the flow of solar plasma much in the fashion 
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suggested by Chapman-Ferraro more than 30 years to explain 
transient variations of terrestrial fields following solar 


activity. 
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A naive and simplified approximation to the very complex 

interaction of the solar plasma impacting the earth's 

magnetic field is shown in Figure 6. Here specular reflection 

from the boundary of the earth's magnetic field is assumed 

for the individual particle trajectories. Within this 

bounding region the dominant magnetic field is that of the 

earth’s and the distortion of the geomagnetic field by currents 

on the containing boundary. This region of space has been 
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referred to as the magnetosphere by Gold since the motion of 
charged particles is completely dominanted by the earth's 
magnetic field. On the assumption of direct impact of the 
solar plasma on the earth's magnetic field, the stand off 
distance of the boundary, depending upon the solar plasma 
momentum flux and geomagnetic field strength, can be predicted. 
Figure 7 summarizes a parameter distribution of stand-off 
distances as a function of particle densities and energies. 

In order to convert from these values to fluxes one multiplies 
particle density by the equivalent proton velocity. The fluxes 
being discussed have been measured by space probes and 
satellites, but this paper does not intend to review this material. 
Interested readers are referenced to the papers by Snyder and 

8 9 10 

Neugebauer, 1963: Bonetti et al . 1963; and Bridge et al . 1964. 

The observed distance of 10.7 R e and estimated plasma velocities 
of 300 to 700 Km/sec leads to plasma densities of 1 to 10 protons 
per cubic centimeter. A general summary of the present field 



environment illustrating the relative position of the radiation 
belts in the magnetosphere is shown in figure 8„ This 
illustration also shows the position in which the Explorer X 
satellite probe mapped the magnetic field of the earth and 
detected only small differences between theoretical and observed 
magnetic fields, j 

The results of the IMP-1 magnetic field experiment for out 

bound orbit number 15 is shown in Figure 9. The similar 

characteristic of a gradual increase in the strength of the 

observed field , over and above that predicted, is clearly seen. 

Beyond this, a region of rapidly fluctuating and low magnetic 

field strength is observed to extend to a distance of 23.1 R e . 

Beyond this distance the magnetic field is reasonably steady and 

stable in orientation. We identify this abrupt temination of 

the fluctuating or tubulent region as a collisionless 

11 

magnetohydrodynamic shock wave associated with the interaction 

of the solar plasma with the earth’s magnetic field. Our 

present understanding of this phenomenon is principally based 

12 

upon an analogy with high speed gas dynamics. In a fluid, 
disturbances can be propagated at a sensible velocity 
characteristic of the medium. In gas dynamics it is the 
acoustic velocity while in magnetohydrodynamics in which ionized 
gases are in motion the appropriate velocity is the Alfven speed. 
This is shown in figure 10 as a function of particle density and 
magnetic field strength. 
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The average magnetic field strength in interplanetary 
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space has been measured to be approximately 5 gammas with extreme 
values between iy and 10 y . Particle densities between 1 
and 10 protons per cm lead to an Alfven velocity in general 
less than 100 Km/sec. Direct measurements of the plasma 
velocities in interplanetary space indicate velocities on the 
order of 400 Km/sec. or more^ and thus the flow of solar plasma 
is super Alfvenic with respect to the propagation of 
characteristic disturbances in the medium. On this basis 
the interaction of the earth's magnetic field with the solar 
plasma becomes much like that of a blunt object in supersonic 
gas flow in which a detaohed shock wave preceeds the body. It 
is separated from the obstructing article, which in this case is 
the magnetosphere boundary and not the earth's surface. A 
summary of the magnetic field environment of the earth in cis- 
lunar space is shown in figure 11. This is based upon the first 
19 orbits of the IMP-1 satellite^ which is the first satellite 
to have performed detailed and accurate measurements of both the 
magnetosphere boundary and the coll ionless shock boundray. 

Three regions of space, 

1 . Magnetosphere 

2. Turbulent Boundary layer and 

3. Interplanetary medium 

are representative of the characteristic aspects of magnetic 
fields which must be simulated on the surface of the ^arth. I The 
I important conclusion deduced from these data is that (very weak 
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magnetic fields are observed in space and must be simulated 

r> 

on the surface of the earth. In addition it also implies 
that the "portable laboratory bench" upon which magnetometer 
sensors are placed must be very clean in the magnetic sense 
so that accurate measurements will be performed. 
geomagnetic Field Representation 

In a region of space containing no sources of magnetic 
flux the magnetic vector F is derivable from a potential V 
which is a solution of Laplace’s equation. In a spherical 
coordinate system (R, 9, (f> ) represent geocentric distance, 

colatitude and longitude east of Greenwich. The potential V 
can be represented in terms of a series of spherical harmonic 
function as 


ao n 


V=R„ 2 2 P„ m (cos0) 

n* I m=0 




1 cos m <p +H l , m sin m <p 


in which 

R e = 6371.2 kilometers, equivalent radius of a spherical Earth 

P^ 111 = the partially normalized associated Legendre polynomials 

introduced by Schmidt and utilized in classical 
14 

geomagnetism. 

G n m , H m = harmonic coefficients for sources internal to 
“ n 

spherical surface (R =R e ) with the dimensions of 


magnetic force 
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G , H m = the harmonic coefficients for sources external 

to the spherical surface (R = R e ) with the 

dimensions of magnetic force. 

From the potential representation the magnetic field F is 

derivable as the negative gradient, F = -grad V, which yields 

straight forwardly the X, Y and Z or southward, eastward and radially 

outward components of the earth's magnetic field. Other 

elements common in geomagnetic field analyses are the horizonal 

and vertical components and declination of the magnetic field. 

There exist at present 10 sets of spherical harmonic 

coefficients for this analytic representation, all differing 
15 - - 

slightly. In general G and H are small (of the order less than 3% of 
G and H) but oecome very important at satellite altitudes 
greater than 5 R e . 

The familiar approximation to the earth's magnetic field 
utilizing a centered dipole is based on the following values 
of harmonic coefficients: 

G-^® = 30470 gammas 

G^ 1 = 3030 gammas and • 45770 

Using only these terms the observed magnetic field can be 
predicted within an accuracy of 87 to 91 per cent at geocentric 
distances of up to several earth radii. These values of 
coefficients correspond to a centered dipole titled 12® with 
respect to the earth's rotational axes and at a longitude 
297.7° east of Greenwich with an equivalent surface magnetic field 
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at the equator of 31,159 gammas. This is obtained from Vestine's 
model A for 1955 which is the best of all simple dipole 
coefficients. It is possible to represent the earth's magnetic 
field by an eccentric dipole which improves upon the accuracy 
and utilizes the first eight coefficients in a spherical 
harmonic expansion. In this representation the centered dipole 
is replaced by a dipole displaced from the center of the earth 
by 436 kilometers and intersecting the surface of the earth 
at 

9 n = 8.1° 9 S = 165.6° 

§ = 275.4° E § = 120. 4°e 

n * s 


Use of this representation yields accuracies of the order of 

90-96% for altitudes less than several R e . At present most 

current students of the geomagnetic field have used coefficient 

sets with less than or equal to 63 elements corresponding to 

n and m being less than 6 (one exception utilized 512 

coefficients) . However, the impact of the solar wind leads 

to a distortion of the earth’s field which at geocentric 

distances greater than 5 contributes a significant portion 

of the observed magnetic field. This has been studied by 
16 

Beard and Mead to determine the contribution of the external 
sources to the magnetic field interior to the magnetosphere. 
They used a centered dipole approximation with an equatorial 
field strength of 31,000 gammas and assumed impact of the solar 
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plasma normal to the dipole axis. Satisfying the boundary 
condition of zero normal component of the net magnetic field 
and assuming a curl free magnetic region interior to the 
magnetosphere leads to the following values for the external 
sources: 

G ° = -0.277 (R e /R b ) 3 gauss equals 27.7 gammas at R fe = 10 R g 

G * = 0.108 (R /Rj 4 gauss equals -1.1 gamma at R. = 10 R 

2 0D u tJ 

At altitudes greater than 5 R g the contribution is large. 

For example at 10 earth radii the external field contributes 

46 gammas while the internal field contributes only 31 gammas. 

At 5 R the contribution of the external field is reduced 
e 

to 37 gammas while the internal field is increased to 248 gammas. 

We see that in the noon meridian the distortion is large, 

approximately 15% at a distance of 5 R g but on the midnight 
meridian the situation is vastly different. Satellites are 
presently investigating this particular region of space and 
theoretical predictions and experimental evidence indicate 
that the earth's magnetic field may trail out far into 
interplanetary space well beyond the orbit of the moon in a 
fashion analogous to the structure of cometary tails in 
interplanetary space. In conclusion, for accuracies of the 
geomagnetic field representation of several percent one must 
restrict altitudes at present to less than about 4 R 0 and include 
external sources for greater altitudes. 
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On the assumption that 


(1) the efflux of solar plasma from the sun is radially 
outward at velocity V 

s 

-6 

(2) the rate of solar rotation is 2.9 x 10 radians 
per second ( =£2 ) 


(3) assuming a photospheric magnetic field strength given 

17 

by B q it is possible, following Parker, to predict 
the strength of the interplanetary magnetic field as 


B=B »(f 

( 2 ) 

where a = radius of sun 

r = heliocentric distance 

This theoretical model of the interplanetary field shows a 
characteristic geometrical configuration in interplanetary 
space which resembles the classical Archimedean spiral in the 
ecliptic plane. A representation of this spiral structure in 
the interplanetary field is shown in Figure 12. The fixed 
rotation rate of the sun and high plasma velocity combine to 
drag out the lines of magnetic force by the highly ionized 
gases and leads to these characteristic structures. In Figure 
13 the theoretical angle, as measured in the plane of the 
ecliptic, which the field makes with a radial line to the sun 
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is shown as a function of plasma velocity. It should be 
noted that the particle density has no bearing on the direction 
of the interplanetary magnetic field strength nor its magnitude. 
It is seen that the theoretical angle <f> is approximately 135° 
near the orbit of the earth. 

Direct measurements of the interplanetary magnetic field 
have been performed from the IMP-1 satellite. These measurements 
indicate a striking consistency with the theoretical model 
just discussed. A sample of the magnetic field measurement in 
interplanetary space obtained on the IMP-1 satellite is shown 
in Figure 14. It is seen that the magnitude is approximately 
5 gammas and that the field is reasonably near the plane of 
the ecliptic and approximately at the azimuthal angle 
theoretically predicted. This indicates, by extrapolation to 
the surface of the sun, that field magnitudes on the order of 
several gauss at least are present. It also means that for 
space probes going inward toward the sun a general increase 
in the magnetic field strength should be observed. Unfortunantly 
the one satellite carrying magnetometers to date which has 

investigated the region of space towards the sun, Mariner II } 

was contaminated by spacecraft magnetic fields and did not allow 

18 

accurate measurements to be performed. 

Simulation of Magnetic Fields in Space 

The preceeding discussion has summarized our present know- 
ledge of magnetic fields in space and presented a foundation 
upon which we can base out requirements for simulation of 
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such magnetic fields. The most common method employed in 
the past has been to utilize circular coils of wire carrying 
currents to create a magnetic field whose magnitude and 
direction can be controlled. The first coil system developed 
was due to Ampere and consisted of a single coil as 
illustrated in Figure 15. in this case the axial component 
of the magnetic field is given by 


H z (0,0,Z)=4- 


2 [(Z-Z 0 ) 2 +A 2 ] 3/2 


(3) 


Two such coil systems placed symmetically with respect to the 
origin at + B give 


H z (0,0,Z)= 2 | [- (z _ B) 4 A 2]3/2 + L(Z+B) 2 +A 2 ] 3/2 


(4) 
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In order to optimize the two coil configuration it is desirable 
to have a maximum field on the Z axis for a minimum current. 
This is obtained analytically by differentiating the above 
formulation, setting it equal to 0 and solving for the various 
values permitted: 


3 H 2 

3z 


= 0 

o,o, z 


(5) 


In addition, in order to provide a maximally homogeneous 
magnetic field one would require a minimum gradient of the 
actual component of the magnetic field. This is determined 
by setting the second derivative of the field with respect to 
Z equal to 0 as shown below 



=0 


0,0, z 


( 6 ) 


Now clearly at Z equals 0 the field is a maximum regardless of 

the value of B as long as the location of the two coils is 

chosen to be symmetric. However, in order that the gradient to 

2 2 

be a minimum then the final result of formula (6) is that 4B = A 
or that A = 2B. This optimum coil configuration utilizing 
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two coils was developed by Helmholtz many years ago and 
has become a standard reference in the simulation of magnetic 
fields in space. 

Maximum uniform magnetic fields utilizing multiple sets 
of coil pairs rely on these two principles: 

1. That pairs of coils symmetrically located with respect 
to the orgin give maximum fields and 

2. Minimum gradients _if appropriately positioned (which 
can be determined by analytical investigations) . 

One can introduce four coils with different numbers of windings on 
each pair of coils or vary the distances and/or size of radii. 

In general the homogeneity of the magnetic field is the most 
important goal since symmetry quite easily achieves the 
maximum field. Indeed a figure of merit can be defined which 
measures the degree of homogeneity as: 


100 


H(x,y,z) 

H(0,0,0)_ 


= % Homogeneity 


(7) 


Now the Helmholtz coil provides a homogeneity of + 1% 
along the axis if the absolute value of Z is less than .32 A. 

It has a maximum error of + 1% off the axis if the absolute value 
of Vx 2 + Y 2 is less than .38 A. It is seen that the region 

of space homogenous in magnitude to 1% is roughly an oblate 
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spheroid whose dimensions of principal axes are given by .32A, 

. 38A , . 38A. 

Now the earth's magnetic field is approximately 50,000 
gammas and one percent of this is 500 gammas so that for 
interplanetary work one needs a uniformity much better than 
one percent. In order to provide a homogeneity of one 

gamma over a working region of space implies a requirement of 
homogeneity of .001%. For a Helmholtz coil this implies a very 
small region of space and a large coil radius. Indeed, for an 
homogeneity of 0.1% the dimensions of the region of space are 
approximately one half the dimensions for 1% accuracy, being 
.16A, .2A and ,2A. For values of .001% the dimensions are 
.04A, .05A and .05A, so that if the radius of the coils is 25 
feet then one has a sphere of approximately 1 foot radius in 
which the field is homogeneous to 1 gamma. Clearly the practical 
liminations of constructing such large coil systems requires a 
re-evaluation of the utility of the very simple Helmholtz 
configuration. It is found that by adding more coils but using 
the same principles governing the development of the Helmholtz 
system one can rapidly improve upon the situation for the 
practical development of coil systems. 

An intriguing possibility but one which has not proven 
successful in practical applications is available in what is 
referred to as a "sine winding" coil system. In this geometry 
the density of turns per unit length actually is held constant 


- *9 ~ 


while the coil is continuously wound on the surface of a sphere. 
This provides a completely uniform magnetic field interior 
to the spherical surface. Unfortunately suc h a configuration 
is not reasonably practical when consideration of access to 
the coil system is made. 

The general approach is to develop hybrid coil systems 

utilizing multiple pairs of coils and employ analytical 

procedures to determine optimum coil constants and configurations. 

One begins by using circular coordinates and expressing the 

potential V in cylindrical harmonics and then determining the 
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appropriate turns and. current values. The important considerations 
in the development of such coil systems are: 

1. The difficulty in implementing the coil manufacture, 
which has lead to consideration in some instances 

of square coil systems, 

2. Current stability and current ratios required, which 
has lead to constant current but separate turns 
ratios for coil pairs rather then attempting to develop 
precise current dividers utilizing resistive ratios 
only; 

3. The accessibility for various geometries, which has lead 
to the development of cubic and square coil systems and 

4. The analytical difficulties in developing optimum coil 
systems, which have lead students of coil geometries 
to consider mainly equal amp terms or sets of coils 
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which lie on the surface of a sphere but at 
different distances - .ong the axis of the system. 

A summary of some existing coil systems and their figure 
of merits is shown in Table I. 


Summary of Some Existing Optimum Coil Systems 


volume of 

Name Dimensions Homogeneity (0.1%) Turns ratio 


Circular Helmholtz 

b = 0.5a 

0.9% 

1:1 

20 

Fanselau 


14% 

equal amp 

21 

Braunbek 


22% 

turns 

Tt M 

22 

McKeehan 


22% 

on surface 

19 

Square Helmholtz 

b = 0.55a 

1% 

of sphere 
1:1 

19 

4 Coil Cubic 

b = ± .33a, +a 

0.8% 

46:43:43:46 

19 

5 Coil Cubic 

b = 0 , ± .5a, ±a 

2.2% 

134:44:66:44:134 


TABLE I 


Presently the Goddard Space Flight Center is fabricating a 
25 foot Braunbek coil system in which three separate sets of 
such coils are employed for three axis cancellation of the 
earth's magnetic field./ The specific parameters of the GSFC 
facility are; 

1. the diameter of the smaller coil is .76388 of the 
larger , 


2 . 


the separation of the smaller is .84565 of the larger 
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diameter, and 

3. the separation of the larger coils is .27803 of the 
larger diameter, and 

4. The real diameter of the coil systems are: 

(a) the vertical is 22 feet, 

(b) the horizonal east-west magnetically oriented 
is 18 feet and 

(c) the horizonal north- south magnetic is 16 feet. 
The orientation of the geomagnetic field at the Goddard 

facility is such that a larger coil system is required for the 
vertical component since it is the largest in magnitude, while 
the remaining coil systems are chosen to be compatible with 
practical implementation of an inter laced set of such coils. 


In Figure 16 is shown a picture of a model of the GSFC coil 
system. The coil constant is 15.4 gammas per milliamp and there are 
36 turns of number 6 copper wire on each coil. The performance 
of the Braunbek coil system is seen, by review of the 
characteristics of Table I, to be equivalent to a McKeehan and 
to have the additional merit that equal amp turns on each coil 
system are required, whereas the other requires different currents 
in each of the coils. The measured performance of the Goddard 
facility provides a homogeneity of .001% in a 1 meter diameter 
sphere with an accuracy of + 0.5 gamma and a resolution of 0.1 


gamma 
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An impbrtant aspect of such coil systems is that they 
be capable of being servo controlled to simulate a stationary 
weak magnetic field in spite of the fact that the earth's 
magnetic field varies as a function of time. A secondary set 
of coils is associated with this coil facility an< i is used 
to monitor the variations of the geomagnetic field. The 
stability of the coil system is approximately 1/2 to 1 gamma and 
the cancellation of the earth's magnetic field can thus be 
automatically controlled. The Goddard facility is not yet 
completely operational although sensible measurements have been 
made utilizing the coil systems. The completed system is scheduled 
for March 1965. As such it will be the largest and most accurate 
facility in the world. 

S pacecraft Applications 

The most important application of any coil system is its 
ability to simulate magnetic fields which represent those to be 
measured in space for calibration of magnetometer sensors. A 
secondary objective is to permit the mapping magnetically of the 
inherent magnetic properties of spacecraft and spacecraft sub- 
systems whose magnetic fields may yield adverse effects on the 
data collected. The following paragraphs summarize briefly 
the specific applications of the knowledge of the earth's 
magnetic field and its interaction with satellites which 
have developed over the past several years. 
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Vanguard I - Spin Decay 

The Vanguard I satellite provided a unique example for 
investigation of the electro magnetic interaction of a spin 
stabilized satellite with the earth's magnetic field over a long 
time scale. The basic physics involved is based upon the fact 
that a rotating sphere in a magnetic field will induce eddy currents 
in a fashion such as to generate currents whose magnetic fields 
not only oppose the spinout whose currents lead to dissipation 
of Kinetic energy. On the assumption that the electro magnetic 
readjustments, as the spin axis changes inclination with respect 
to the magnetic field, occur at a rate much greater than the 
mechanical readjustments, the theoretical development and expla- 
nation 0 f fhe Vanguard I spin rate decay can be understood. 

The first investigation was conducted by j. p. Vinti assuming 
a spherical shell of uniform conductivity and non-magnetic 
material* The characteristic summary of the results is that 
the induced eddy currents create a torque on the spinning 
satellite proportional to the inertial moment of the satellite, 
and the angular spin rate and the magnetic field strength squared 
as shown in the formulae below: 


Torque a I(wxB)xB 

( 8 ) 


where I is the angular moment of inertia 
w is the angular velocity 

B is the magnetic field vector 
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The net result of the torque can be separated into a spin 
decay which j.s proportional to the angular spin rate and its 
product with the magnetic field perpendicular to the 
spin rate as shown in the following formula: 

L (decay) a wB ± 2 

Similarly the precession torque is 


L (precession) a w Bj, B 1 1 

The Vanguard I satellite was launched March 17, 1958 into 

an orbit with elevations between 650 to 3,900 Km at which the 

earth's magnetic field was approximately 40 to 30 thousand 

gammas. The spin decay studies indicated the satellite slowed 

24 

down from an initial rate of 3.0 rps to 0.05 rps in 2.5 years. 
These measurements were performed by investigation of the periodic 
fading of radio signals associated with the antenna pattern of 
the satellite. The first anomaly was detected in December 1958 
and a rapid evaluation of theoretical models which go back to 
Hans Hertz in 1896 indicated that classical EM theory could 
explain the observed anomaly. 

Injun HI Passive Magnetic Orientation 

The motions of charged particles trapped in the earth's 
magnetic field are generally helical spirals about field lines 

with associated drift of electrons east and protons west due 
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to the spatial gradient of the earth’s magnetic field. The 
particles "bounce" back and forth along 

a line of force between their "mirror" points. The angle that 

the spiral makes with the field lines is called the pitch angle and 

and is an important parameter to be measured in charged 

particle investigations in space. Since it is known that the 

first adiabatic invariant of particle motion is preserved then 

measuring the pitch angle distribution of the particles at a 

point in space predicts the particle flux characteristics along 

that line of force. The variation in pitch angle is shown by the 

following formula in which iB Q is the field strength when a - 0° 

and B is the field strength at the mirror point when a = 90°. 
m 


B 0 = B m sin 2 a 90 o 


and in general 


sin 2 a = 



( 10 ) 


In order to separate by direct detector analysis the pitch 
angle distribution requires unique instrumentation. An ingeneous 

suggestion has resolved this difficulty by placing a bar magnetic 
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in the spin stabilizated spacecraft Injun III so that the 
satellite spin axis becomes parallel to the local magnetic field. 
Recalling the results of the previous discussion on Vanguard I, 
it is found that both the decay and the precession torques are 
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zero for perfect alignment of the magnetic vector and spin 
vector of the satellite. In the Injun III satellite an 
Alnico-V magnet 22 inches long and 1 inch square was employed 
and permalloy rods perpendicular to the spin axes were 
employed to damp the periodic motion induced by the variable 
aspect to the earth’s magnetic field. The Injun III 
satellite was launched December 13, 1962 and achieved an 
apogee of 2787 Km. and a perigee of 237 Km. The natural 
period of this satellite was approximately two minutes 
and the satellite damped to being almost perfectly aligned 
with the earth's magnetic field in approximately three days. 

The success of this particular passive orientation device 
was checked by onboard magnetometers measuring the magnetic 
field perpendicular to the spin axis of the satellite. 

Tiros-Active Attitude Control System 

The Tiros wheel satellite, number 9, will be launched into 
a circular orbit of 13,000 Km with a unique active attitude 
control system utilizing the principles already discussed. It 
will be placed into a near polar orbit with the spin axis 
perpendicular to the orbital plane and sun synchronized so 
that the precession of the orbital plane due to the earth's 
equatorial bulge will maintain the satellite in a solar 
orientation so that a constant illumination of the observed 
earth's surface is maintained. The magnetic aspect control 
system (MASC) will keep the spin within its proper limits and 
the quarter orbit magnetic attitude control system (QOMAC) will keep 
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the spin axis perpendicular to the orbital plane. These 
systems are illustrated in figure 17 with respect to orientation 
to the satellite spin axis 

The need for such current carrying coil systems was not 
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recognized until after Tiros I. On successive Tiros, numbers 
II through VIII, the magnetic attitude control system was 
utilized mainly to provide proper cancellation of the satellites 
magnetic field so that spin torques on the satellite were 
negligible. It is possible to alter the polarity of the space- 
craft dipole moment by these coil systems and indeed the 
magnetic properties of the satellite are planned to be cancelled 
continously by these attitude systems. 

SUMMARY 

The above discussion hhs been very brief but illustrates 
the utilization of the earth's magnetic field in attitude 
control systems going from strictly interaction to highly 
sophisticated computer programmed active attitude systems. The 
necessity for verifying the operation of such systems is clear 
and the utilization of large coil facilities for attitude studies 
such as that being fabricated at GSFC is clearly indicated. 
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LIST OF FIGURES 


1. Tabular summary of US and USSR earth satellites launched to 
date which have provided measurements of the geomagnetic 
field. Included are a brief description of the type of 
instrumentation employed and the associated dynamic range 
as measured in gammas. The sensitivity is normally 
determined as a percentage of the maximum range (frequently 1%) 
except in the cases of those satellites for which large 
spacecraft magnetic fields were present such as on the 
Russian spacecraft Lunik I and II. The distance over which 

sensible magnetic field measurements were performed is 
indicated in units of earth radii, (see Figure 2) 

2. Tabular summary of earth satellites launched to date measuring 
the geomagnetic field. Indicated are the launch dates, 

the inclination of the orbital plane to the earth's equator 
(when this parameter is significant) and the life time of 
the satellite as measured inchys. Note that some spacecraft 
clearly traversed the indicated region only once. 

3. Tabular summary of US space probes launched to date which 
have provided measurements of the interplanetary magnetic 
field. Prior to the launch of the IMP-I spacecraft no 
accurate and precise measurements of the interplanetary field 
had been performed due to various limitations as indicated in 
the tabular summary under the heading "comments". 



2 


4. Results of the Explorer X measurements of the geomagnetic 
field from 4 to 12.5 R 0 . The observations are indicated by 

solid dots connected by straight line segments, the 
theoretical values indicated by solid lines and interpreted 
fit assuming a constant field indicated by the dashed line. 

5. Results of the traversal of the magnetosphere boundary by 
IMP- I on inbound orbit pass No. 1. Observed values are 
indicated by open and closed circles connected by straight 
line segments as provided by both the rubidium vapor magne- 
tometer and the fluxgate magnetometers carried onboard the 
satellite. The dashed curves represent the theoretical 
magnetic field extrapolated from terrestrial surface 
measurements. At these distances the only contributing 
factor in the multipde expansion of the earth’s field is 
the dipole moment of the earth. The abrupt transition at 
10.7 R_ identified as the boundary of the magnetosphere, the 
magnetopause . 

6. Naive representation of the interaction of the solar plasma 
with the geomagnetic field. Direct impact of the plasma 
with the magnetic field is shown as being specularly 
reflected from the geomagnetic boundary. The distance to 
the boundary at the subsolar point is given by: 


R 

m 


R 

e 



4rr mn V ^ 
s 
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where R i s the radius of the earth, B is the equatorial 

u 

magnetic field strength and V the velocity of the solar 

b 

3 

plasma with density n p/cm . 

7. Theoretical size of the magnetosphere at the subsolar point 
assuming normal impact of the solar plasma on the geomagnetic 
field. For a magnetosphere of size 10.7 R at the 
stagnation point, and assumed velocity of 400 km/sec 

the deduced equivalent plasma density is 2 protons per cubic 
centimeter. 

8. Summary pre- IMP-1 illustration of the confined geomagnetic 
field of the earth and the location of the radiation belts. 

The flow of solar plasma is taken to be directly from the 
sun and the trace of the Explorer X trajectory is shown 

on the magnetic meridian plane. 

9. The results of the magnetic field experiment carried on the 
IMP- I spacecraft from orbit No. 15, January 21, 1964. These 
data illustrate the outbound traversal of the magnetosphere 
boundary at 15.7 R g and the collisionless magnetohydrodynamic 
shock wave at 22.7 R e , with a possible precurser occurring 

at 24.5 R . 

e 

10. Alfven magnetohydrodynamic phase velocity of wave propagation 
as a function of magnetic field strength and plasma density. 
Representative values for the interplanetary medium are chosen 
in this diagram. 
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11. Comparison of the observed positions of the boundary of the 
magnetosphere and shock wave as shown by solid dots with 
the theoretical positions according to the computations of 
Spreiter and Jones (1963) adjusted for a different gas 
dynamic specific heat capacity ratio. Very good agreement 
is obtained with this modification of their treatment and an 
aberration by 5° to accommodate the heliocentric orbital motion 
of the earth around the sun. The distance to the magnetosphere 
boundary at the subsolar point is 10.25 R e and the distance 

to the shock wave boundary is 13.4 R^. 

12. Schematic illustration of the spiral interplanetary magnetic 
field as obtained by assuming a uniformly expanding, 
infinitely conducting solar corona. The topology of the 
interplanetary magnetic field changes significantly as the 
velocity of the plasma increases. At high velocities the 
field lines are radially directed, paralleling the flow field 
from the sun, while for very low velocities the field lines 
are tightly wound in a Archemedian spiral configuration 

as shown on the left hand side. 

13. Theoretical direction of the interplanetary magnetic field 
as measured by the azimuthal angle 0 in the plane of the 
ecliptic as a function of the solar wind velocity V . 

14. Representative measurements of the interplanetary magnetic 
field from IMP-I during orbit 15, January 22, 1964. These 
data follow immediately after those of Figure 9 and indicate 
the characteristic of the interplanetary magnetic field 
magnitude and variability of the direction. 
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15. Geometry and relative scale of the Ampere coil and Helmholtz 
coil configurations utilized to nullify the earth's magnetic 
field in space simulation. The Z axis of the coil system 

is normal to the plane of the coils. 

16. Model of the Goddard Space Flight Center, Magnetic Fields 
component test facility located adjacent to the laboratory 
in Greenbelt, Md. This large scale coil facility will be 
utilized for calibration of magnetometer sensors and mapping 
of magnetic field properties of spacecraft and subsystems. 

17. Schematic illustration of the Tiros wheel satellite 

geometical configuration of coil systems employed to control 

the aspect and attitude of the satellite. Electrical currents 

can flow in either polarity on the coils as indicated under 
ground computer control* 
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